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INTRODUCTION 


Optical  imaging  through  a  scattering  atmosphere  such  as  fog  and 
smoke  has  been  an  unresolved  problem  since  ancient  times.  The 
compass  was  invented  by  the  Chinese  centuries  ago  to  allow  their 
army  to  advance  in  heavy  fog.  The  interest  in  the  problem  of 
atmospheric  optical  transmission  is  heightened  by  the  development  of 
optical  communication,  LADAR  and  laser  weapon  technologies.  A  great 
deal  of  effort  has  been  geared  toward  the  phase  distorting  aspect  of 
atmospheric  turbulence.  However,  there  have  been  few  viable 
solutions  to  the  problem  of  scattering. 

Attempts  to  recreate  the  target  image  with  the  scattered  light 
are  generally  futile  since  it  is  almost  an  impossible  task  to  deter¬ 
mine  the  time-varying  impulse  response  of  a  scattering  medium.  One 
may  however,  try  to  separate  the  components  contributed  by  the  scat¬ 
tered  and  unscattered  light  fields  and  then  reconstruct  the  target 
image  with  the  unscattered  light.  Holography  has  been  proposed  at 
various  times  as  a  means  of  separating  the  scattered  and  unscattered 
light  fields.  The  effectiveness  of  holographic  techniques  in  imaging 
through  fog  has  been  demonstrated  quite  successfully  in  the  labora¬ 
tory.  However,  the  conventional  holographic  techniques  that  were 
employed  are  extremely  sensitive  to  vibration  and  atmospheric  turbu¬ 
lence.  This  low  tolerance  renders  them  impractical  in  realistic 
field  conditions.  In  this  report,  we  shall  introduce  a  system  that 
utilizes  an  imaging  grating  interferometer  to  produce  a  local  refer¬ 
ence  beam  hologram  with  the  unscattered  light  waves.  Unlike  conven¬ 
tional  holographic  arrangements,  this  proposed  concept  does  not  re¬ 
quire  the  use  of  a  separate  reference  beam.  This  enables  the  con¬ 
struction  of  a  small  and  compact  system.  More  importantly,  besides 
providing  significant  improvement  in  imaging  through  a  highly  scat¬ 
tering  atmosphere,  the  proposed  system  possesses  much  higher  toler¬ 
ance  to  vibration  and  atmospheric  turbulence,  distinguishing  it  as  a 
system  feasible  for  field  applications. 


In  the  next  section,  we  define  the  general  problem  and  establish 
the  rationale  for  utilizing  holographic  techniques  in  imaging 
through  scattering  media.  The  use  of  the  conventional  holographic 
technique  is  described  in  Section  3  and  its  drawbacks  in  field 
applications  are  explained.  In  Section  4,  the  grating 
interferometric  imaging  system  is  introduced.  Its  use  in  imaging 
through  scattering  media  is  analyzed  and  its  performance  is 
estimated.  The  optimal  design  and  fabrication  technique  for  the 
grating  system  is  given.  Experimental  results  with  the  grating 
interferometric  system  for  imaging  through  scattering  media  are 
presented . 
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BACKGROUND 


When  the  light  field  reflected  by  the  target  propagates  through 
the  scattering  atmosphere,  it  is  attenuated  by  absorption  and  scat¬ 
tering.  In  addition,  the  scattered  light  creates  a  time  varying 
background  light  field.  Often,  the  intensity  of  this  background 
light  is  substantially  higher  than  the  attenuated  signal  light.  The 
result  is  an  image  of  drastically  reduced  contrast.  If  the  image 
contrast  is  below  2%,  then  the  target  image  would  not  be  recogniz¬ 
able  by  a  human  observer^.  Contrast  is  defined  as 

x  ioo» 

xmax  Amin 

where  I  and  I  ■  are  the  maximum  and  minimum  intensity  levels 
max  min 

of  the  image  received  by  the  observer.  If  the  target  has  an  intrin¬ 
sic  contrast  of  100%,  then  the  intensity  of  the  unscattered  light 
would  correspond  to 

*  unseat  =  *max  -  Jmin 


Tnus,  the  target  image  would  not  be  recognizable  by  a  human  observer 
if  the  SNR  of  the  image  is  lower  than  1:25. 

The  farther  the  target  is  from  the  observer,  the  larger  the 

volume  of  scatterers  that  the  reflected  light  must  propagate 

through.  The  result  is  an  exponential  decrease  in  image  contrast 

with  target  range.  The  image  contrast  C,  as  a  function  of  target 

2  1 

range  R  can  be  expressed  as 

Cj  =  CQ  e~oR  (2-1) 

where  CQ  is  the  target  contrast  observed  at  close  range  and  a  is 
the  visibility  attenuation  coefficient.  The  coefficient  is  dependent 
on  the  wavelength  of  the  illuminating  light  source.  The  amount  of 
scattering  tends  to  decrease  with  longer  wavelengths. 

The  average  visibility  attenuation  coefficient  for  moderately 
heavy  fog  under  sunlight  illumination  for  example,  is  about  10 
km-1.  Thus  for  a  range  of  0.35  km,  the  image  contrast  of  a  target 
with  100  %  intrinsic  contrast  would  be 

Cj  =  100  x  10~3'5  =  3% 

The  target  image  would  be  barely  recognizable  by  a  human  observer. 
However,  the  fact  that  a  target  is  located  within  the  visibility 
range  does  not  necessarily  imply  that  it  will  be  visible.  The  2% 
minimum  contrast  criterion  applies  only  to  a  simple  two-tone  target. 
For  a  complicated  continuous-tone  target  that  is  imbedded  in  a  back¬ 
ground,  the  image  contrast  required  for  human  recognition  would  be 
substantial ly  higher. 

In  order  to  separate  the  components  contributed  by  the  scattered 
and  unscattered  light,  we  must  look  for  a  physical  characteristic 
where  the  two  light  fields  differ.  The  problem  is  much  simplified 
if  we  can  make  the  assumption  that  the  target  is  stationary  within 
the  detection  time  and  the  scatterers  are  in  constant  motion.  The 
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scatterers  are  either  moving  in  bulk  as  they  are  carried  by  air  cur¬ 
rent  or  moving  in  a  random  manner  due  to  Brownian  motion.  First,  we 
examine  the  difference  in  velocities  between  the  target  and  the 
scatterers.  If  the  target  is  illuminated  with  monochromatic  light, 
the  light  reflected  from  the  stationary  target  is  doppler  shifted 
when  it  is  scattered  by  the  moving  scatterers.  The  bulk  type  motion 
would  result  in  a  shift  in  the  center  frequency  of  the  light  and  the 
random  motion  will  cause  a  broadening  of  the  spectral  bandwidth. 
With  the  use  of  a  narrow  band  spectral  filter  such  as  a  Fabry  Perot 
interferometer,  it  would  be  possible  to  separate  out  the  unscattered 
light  from  the  doppler  shifted  scattered  light.  However,  this  scheme 
has  some  drawbacks.  High  wind  speed  can  not  be  relied  on  and  if  the 
wind  is  blowing  traversely  across  the  light  path,  no  doppler  shift 
will  be  observed.  Brownian  motion  is  always  present  but  unless  the 
scattering  particles  are  very  small,  the  particle  velocities  would 
not  be  adequate  to  produce  significant  doppler  shift.  Moreover,  even 
though  the  particles  are  in  constant  motion,  their  mean  velocity  is 
zero.  Thus  a  substantial  amount  of  scattered  light  would  remain 
within  the  spectral  band  of  the  light  source  and  the  separation  be¬ 
tween  the  two  light  fields  cannot  be  complete.  Spectral  filtering 
with  a  Fabry  Perot  etalon  also  would  not  produce  an  image  directly 
since  the  system  only  provides  the  information  on  the  power  of  the 
incoming  light  field  that  is  within  its  passband.  To  produce  an 
image,  it  would  be  necessary  to  look  at  a  small  area  at  a  time  and 
build  up  into  a  complete  image  by  some  form  of  scanning  mechanism. 

Alternatively,  one  can  make  use  of  the  fact  that  the  moving 

scatterers  create  a  time-varying  light  field  for  the  scattered 

3-4 

light.  With  an  achromatic  grating  interferometer  for  example, 
the  lateral  position  of  a  point  radiator  is  mapped  into  the  spatial 
frequency  of  a  sinusoidal  fringe  pattern.  The  moving  scatterers 
would  therefore  create  fringe  patterns  of  varying  frequencies.  If 
the  moving  scatterers  are  constantly  moving  to  new,  never  occupied 
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positions  in  the  field  of  view  during  the  integration  time  T,  the 
time-averaged  contrast  of  the  fringes  formed  by  the  scattered  light 
would  only  be  1/T  of  that  produced  by  the  stationary  point  target. 

However,  this  would  approximately  be  true  only  if  the  number  of 
scatterers  is  much  smaller  than  the  space  bandwidth  in  the  field  of 
view  of  the  Imaging  system.  In  most  practical  applications  (e.g., 
imaging  through  fog),  the  sizes  of  the  scatterers  and  the  spaclngs 
between  them  are  both  beyond  the  resolution  limit  of  the  imaging 
system.  For  an  incoherently  illuminated  target,  the  modulus-squared 
of  the  point  spreads  functions  will  merge  to  form  a  continuous  inten¬ 
sity  pattern.  Thus,  even  though  the  scatterers  are  in  constant  mo¬ 
tion,  the  intensity  distribution  of  the  scattered  light  field  would 
appear  stationary  to  the  Imaging  system.  The  SNR  cannot  be  improved 
by  time-averaging  if  the  target  is  incoherently  illuminated. 

In  order  to  separate  the  scattered  and  unscattered  light,  it 

would  be  necessary  to  illuminate  the  target  by  coherent  light.  If 

the  target  is  imaged  through  a  scattering  media,  the  unscattered 

5 

light  would  form  a  speckled  image  . 

|f(x,  y)  *  h(x,  y)|2 

where  f(x,  y)  is  the  target  function  and  h(x,  y)  is  the  impulse  re¬ 
sponse  of  the  imaging  system  and  *  denotes  convolution  operation. 

On  the  other  hand,  the  scattered  light  produces  a  time-varying 
speckle  pattern  which  in  time  averages  out  into  a  uniform  bias. 
Thus,  by  recording  the  image  over  a  sufficiently  long  integration 
time,  we  would  obtain 

Ta(x,  y)  -  B2  ♦  A2 |f (x,  y)  *  h(x,  y)|  2  (2-2) 

2  2 

where  B  and  A  are  the  intensities  of  the  scattered  and  unscat¬ 
tered  light.  Inserting  the  transparency  with  an  amplitude  trans¬ 
mittance  of  Ta(x,  y)  Into  the  optical  processor  as  shown  in  Figure 
1,  we  would  have,  at  the  spatial  frequency  plane5. 
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9f[TA(x,  y)]  ,  b2«  +  A2[F(p,  q)  H(p,  q)  *  F(-p,  -q)  ff(-p,  -q)](2-3) 

where  the  top  bars  denote  conjugates  and  <0T(  • )  is  the  Fourier 
transform, 

00  00 

_  ff  -i(px+qy) 

3T[Ta(x,  y)]  =  JJ  Ta(x,  y)  e  dx  dy  (2-4) 

•00.00 

For  a  diffuse  object  |  F(p,  q)|  =  constant  and  the  spectral  distri¬ 
bution  of  the  target  Image  would  simply  equal  to  |H(p,q)*H(p,q) |  as  il¬ 
lustrated  in  Figure  2.  Since  the  spectrum  of  the  uniform  bias  con¬ 
centrates  at  the  origin  of  the  frequency  plane,  it  can  be  filtered 
out  with  a  dc  stop.  At  the  output,  we  would  obtain  a  target  image 
of  improved  contrast.  However  the  improvement  achieved  with  this 
simple  technique  is  limited.  First  of  all,  the  Fourier  transform  is 
performed  over  a  finite  aperture.  Thus,  instead  of  a  delta  function, 

we  would  obtain  for  the  bias  term  a  diffraction  pattern  corresponding 

2  2 

to  the  Fourier  transform  of  the  input  aperture.  Since  B  »  A  , 
the  sidelobes  of  the  dc  term  could  contribute  signif icantly  to  the 
filter  image.  Secondly,  the  energy  of  the  image  spectrum  also  cen¬ 
ters  around  the  origin.  In  performing  the  high  pass  filtering,  a 
substantial  amount  of  signal  energy  would  be  lost.  To  make  the  mat¬ 
ter  worse,  system  noise  such  as  film  grain  noise  is  highest  at  low 
frequencies®.  The  high  noise  level  together  with  the  inefficient 
filtering  process  severely  limit  the  amount  of  contrast  improvement 
attainable. 

To  produce  a  much  higher  degree  of  contrast  improvement,  an 
effective  approach  is  to  shift  the  signal  image  spectrum  away  from 
the  dc  term  to  a  region  where  the  system  noise  is  lower.  This  would 
also  allow  the  entire  image  spectrum  to  pass  through  the  filter.  To 
produce  the  necessary  carrier  frequency,  various  holographic  techni¬ 
ques  can  be  utilized. 
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FIGURE  1.  OPTICAL  PROCESSING  OF  COHERENTLY  RECORDED  TARGET  IMAGE 


FIGURE  2.  SPECTRAL  DISTRIBUTION  OF  RECORDED  IMAGE 
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CONVENTIONAL  HOLOGRAPHIC  IMAGING  TECHNIQUES 


If  the  Incoming  light  field  Is  made  to  Interfere  with  a  station¬ 
ary  coherent  light  beam  (reference  beam),  the  unscattered  portion  of 
the  received  light  field  will  produce  a  stationary  fringe  pattern 
while  the  scattered  light  will  produce  time-varying  patterns.  Re¬ 
cording  these  Interference  patterns  with  an  Integrating  detector  such 
a  photographic  film  or  CCO  detector  arrays,  the  time-varying  patterns 
will  in  time  be  averaged  Into  a  uniform  bias.  Except  for  this  added 
background  bias,  the  time-averaged  Interference  pattern  would  be 
identical  to  that  produced  by  the  object  light  field  without  the 
presence  of  the  scattering  media. 


With  the  addition  of  a  reference  beam  and  the  recording  of  the 
interference  pattern  as  illustrated  In  Figure  3,  a  hologram  of  the 
target  is  being  recorded  through  the  scattering  medium/  The  light 
intensity  distribution  at  the  recording  plane  can  be  written  as5 


I(x,  y)  -  B2  ♦ 


.  IkR 
e 

ITT 


-e)2] 


da  dB 


-Ik  sin  ex 


♦  C  e 


(3-1) 


where  f(a,  e)  Is  the  light  distribution  at  the  target  plane  and  R  is 
the  distance  between  the  target  and  the  recording  plane. 

After  recording  this  light  field  linearly,  the  developed  hologram 
is  then  illuminated  by  an  off  axis  coherent  light  beam.  The  field 
of  the  unscattered  light  will  be  reconstructed  to  product  the  light 
field  of  the  unscattered  light.  By  Imaging  this  light  field  (or  by 
simply  looking  Into  the  reconstructed  beam),  the  target  image  can  be 
obtained  with  much  of  the  background  bias  eliminated.  However,  with 
the  cases  of  interest,  ^  ^  must  maximize  the  signal 

level  to  produce  a  satisfactory  Image.  Since  the  point  radiators  on 
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FIGURE  3.  RECORDING  AND  RECONSTRUCTION  OF  H0L06RAPI 

A  SCATTERING  MEDIUM 
10 


the  Illuminated  target  are  all  concentrated  at  a  lateral  plane  while 
the  scatterers  are  distributed  over  a  large  volume,  the  signal  level 
can  be  maximized  by  Imaging  the  target  on  the  recording  plane  and 
construct  an  Image  plane  hologram  as  Illustrated  In  Figure  4.  The 
light  intensity  distribution  at  the  recording  can  then  be  written  as 


A  |f(x,  y)  *  h(x,  y) |  e 


1*(x,y) 


-  Ik  sin  ex 


+  Ce 


(3-2) 


where  A  and  B  are  the  amplitudes  of  the  unscattered  and  scattered 

light,  C  is  the  amplitude  of  the  reference  beam,  f(x,  y)  *  h(x,  y) 

is  the  amplitude  and  #(x,  y)  is  the  phase  distributions  of  the  image 

field.  For  maximum  diffraction  efficiency,  the  amplitude  of  the 

reference  beam  should  be  made  to  be  the  same  as  the  unscattered 

2  2 

light.  That  is,  A  *  C.  And  since  B  »  A  for  the  cases  of 
interest,  the  light  intensity  distribution  at  the  image  plane  can  be 
approximated  by 


I(x,  y)  =  1  ♦ 


f(x,  y)  *  h(x,  y) 


•cos  [k  sin  ex  +  #(x,  y)]  (3-3) 

Recording  the  intensity  pattern  linearly  on  a  photographic  plate  and 
then  illuminating  the  developed  transparency  with  the  reference  beam, 
the  unscattered  object  wave  would  be  reconstructed.  The  recon¬ 
structed  wavefront  can  then  be  imaged  onto  the  output  plane  with  a 
relay  lens  as  shown  In  Figure  5.  We  would  once  again  obtain  the  tar¬ 
get  image  with  a  substantially  improved  contrast. 

Such  a  technique  has  been  demonstrated  in  the  laboratory  success¬ 
fully.  However,  its  effectiveness  is  severely  limited  when  applied 
in  actual  field  conditions.  First  of  all,  to  image  a  target  a  dis¬ 
tance  away,  an  equivalent  path  length  must  be  provided  for  the 
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FIGURE  4.  RECORDING  OF  A  IMAGE  PLANE  HOLOGRAM 


FIGURE  5.  OPTICAL  PROCESSING  OF  IMAGE  PLANE  HOLOGRAM 


reference  beam.  This  would  require  the  prior  knowledge  of  the  target 
range  to  within  the  coherence  length  of  the  light  source.  Since  the 
target  is  obscured,  such  a  precise  knowledge  of  the  target  range 
generally  cannot  be  assumed.  Secondly,  the  ideal  optical  arrangement 
depicted  in  Figure  4,  is  not  a  realistic  simulation  of  what  could  be 
implemented  in  the  field.  For  example,  if  the  light  source  is  lo¬ 
cated  near  the  receiver  as  shown  in  Figure  6,  there  would  be  severe 
backscattering  from  the  illuminating  bean.  Even  if  we  assume  that 
this  backscattering  can  be  eliminated  by  pulse  gating  techniques, 
the  illuminating  beam  still  has  to  propagate  through  the  scattering 
atmosphere  to  the  target.  Only  the  portion  of  the  illuminating  light 
that  travels  to  and  back  from  the  target  unscattered  will  be  able  to 
produce  stationary  fringes.  The  light  intensity  distribution  ob¬ 
tained  with  such  an  arrangement  can  be  approximated  as 

2 

I(x,  y)  -  1  ♦  2^  |f(x,  y)  *  h(x,  y)| 

cos  [k  sin  ox  +  0(x,  y)]  (3-4) 

On  the  other  hand,  if  the  light  source  is  placed  near  the  target 
as  illustrated  in  Figure  7,  the  reference  beam  would  be  scattered  as 
it  travels  to  the  receiver,  creating  additional  forward  scattered 

light.  Under  optimum  beam  ratio,  the  light  intensity  distribution 
would  be 

I(x,  y)  -  1  |f(x,  y)  *  h(x,  y)| 

cos  [K  sin  ox  ♦  $(x,  y)]  (3-5) 

It  can  be  shown  that  the  SNR  of  the  reconstructed  image  is  pro¬ 

portional  to  where  M  Is  the  modulation  depth  of  the  recording. 
Assuming  that  the  recording  is  linear,  then  the  ratio  of  the  SNR 


between  the  ideal  condition  of  Figure  4  and  the  more  realistic  con¬ 
ditions  of  Figure  6  where  the  light  source  is  placed  near  the  re¬ 
ceiver  and  Figure  7  where  the  light  source  is  placed  near  the  target 
would  be 


2  2 

For  example,  if  A  /B  =  1/1 00 ,  then  the  ratio  would  be 

i  1  .  1 

•  T0Cf  *  ? 

Thus,  a  much  better  SNR  can  be  obtained  with  the  light  source  placed 
near  the  target  instead  of  the  receiver.  However,  the  SNR  would 
still  be  a  factor  of  4  lower  than  the  ideal  condition  where  the  light 
beams  can  illuminate  both  the  target  and  the  receiver  without  first 
being  scattered.  It  should  be  noted  here  that  the  ideal  condition 
can  be  achieved  by  propagating  a  narrow  reference  beam  to  the  re¬ 
ceiver  and  then  expanding  the  beam  near  the  receiver  with  a  spatial 
filter  as  shown  in  Figure  8.  The  narrow  beam  will  produce  very  lit¬ 
tle  additional  scattered  light  and  the  spatial  filter  will  allow  only 
the  unscattered  light  to  pass  through.  However,  the  receiver  would 
also  be  obscured  by  the  scattering  atmosphere  and  it  would  be  diffi¬ 
cult  to  aim  a  narrow  beam  precisely  at  the  receiving  optics  in  actual 
field  conditions.  Such  an  arrangement  would  be  restricted  to  appli¬ 
cations  where  the  optics  can  be  fixed  in  prealigned  positions. 

The  most  important  drawback  of  the  conventional  holographic  tech¬ 
nique  however,  is  its  extreme  sensitivity  to  phase  shift  between  the 
interfering  beams.  Such  a  phase  shift  could  be  due  to  displacements 
between  the  various  components  or  changes  in  the  refractive  index  of 
the  medium  in  the  light  path.  In  a  laboratory,  all  the  components 
including  the  target,  receiving  optics  and  light  source  can  be 
mounted  rigidly  on  an  isolated  optical  bench.  In  a  controlled  labo¬ 
ratory  environment,  the  refractive  index  of  the  air  within  the  light 
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paths  can  also  be  maintained  at  a  very  constant  level.  However,  in 
real  world  application  of  the  holographic  system,  the  various  compo¬ 
nents  controlling  the  interfering  beams  would  be  spaced  widely  apart 
with  a  dynamic  air  mass  between  them.  These  components  would  vibrate 
independently  and  it  is  likely  that  the  relative  positions  of  the 
components  would  change  more  than  a  wavelength  during  the  integration 
time.  Moreover,  with  the  interfering  light  beams  traveling  through 
different  light  paths,  atmospheric  effect  such  as  turbulence  could 
easily  change  the  phase  difference  between  the  beams  by  a  wavelength 
or  more.  Thus  even  though  the  conventional  holographic  technique 
has  been  applied  successfully  in  the  laboratory,  it  cannot  be  imple¬ 
mented  effectively  outside  the  laboratory  under  actual  field  condi¬ 
tions.  The  stability  requirement  for  the  conventional  holographic 
technique  will  be  discussed  in  better  details  in  a  later  section. 
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4 

LOCAL  REFERENCE  BEAM  HOLOGRAPHIC  TECHNIQUE 


The  basic  concept  employed  by  the  holographic  technique  is  to 
create  sinusoidal  modulation  for  unscattered  light.  The  modulated 
signal  is  then  separated  from  the  dc  bias  created  by  the  scattered 
light  using  high  pass  or  band  pass  spatial  filtering.  The  need  for 
a  reference  beam  however,  has  limited  their  usefulness  in  practical 
applications.  As  we  shall  show,  the  same  basic  concept  can  be  imple¬ 
mented  without  the  use  of  a  separate  reference  beam.  We  can  produce 
the  sinusoidal  modulation  of  the  unscattered  light  by  interfering 
the  image  light  field  with  itself. 

4.1  DUAL-APERTURE  SPECKLE  INTERFEROMETER 

Let  us  first  consider  the  simple  system  illustrated  in  Figure 
9.  It  consists  of  an  imaging  lens  with  two  apertures  of  diameter  d 
and  separation  D.  For  a  target  f(x,  y),  we  would  have  at  the  image 
plane,  a  light  distribution  of 


where  h(x,  y)  is  the  impulse  response  of  the  imaging  system. 

We  have  therefore  obtained  a  sinusoidally  modulated  target  image 
whose  spatial  frequency  is  proportional  to  the  separation  distance  D 

O 

between  the  apertures  .  If  a  scattering  medium  is  placed  between 
the  target  and  imaging  lens,  the  light  scattered  by  the  medium  will 
lose  its  spatial  coherence  and  the  scattered  light  fields  will  add 
by  their  intensities  while  the  unscattered  light  fields  interfere 
with  each  other.  The  resulting  intensity  distribution  can  be  written 
as 

B2  +  A2  | f ( x,  y)  *  h(x,  y)|2  jl  +  cos^(M^)  +  *(x,  y)jj  (4-2) 


I(X,  y)  =  | f ( X,  y)  *  h(x,  y)| 


(1  +  cos 


2irDx 

X 


>]! 


«*.  y)  (4-D 
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2  2 

where  B  and  A  are  the  intensity  of  the  scattered  and  unscat¬ 
tered  light  at  the  image  plane  respectively.  Thus,  aside  from  a  de¬ 
crease  in  fringe  visibility,  the  light  pattern  is  essentially  un¬ 
changed  by  the  scattering.  If  this  light  pattern  is  recorded  linear¬ 
ly  on  a  photographic  film,  the  modulated  target  image  can  be  ex¬ 
tracted  from  the  background  bias  by  simple  bandpass  filtering  using 
an  optical  processor.  To  optimize  SNR,  the  filter  bandwidth  should 
be  equivalent  to  the  signal  bandwidth.  The  bandwidth  of  the  modu¬ 
lated  term  is  equal  to  the  width  of 

_  2 

<7lf(x.  y)  *  h(x,  y)  |  =|  F(p,  q)  H(p,  q) 

*  F(-p,  -q)  H(-p,  -q) | 

Once  again  for  a  diffused  Target  we  may  assume  that  |  F(p,  q) |  -  con¬ 
stant  and  the  image  bandwidth  would  be  determined  by 

|H(p,  q)  *  H(-p,  -q) | 

Thus,  the  image  bandwidth  would  be  equal  to  2d/\i,  the  width  of  the 
incoherent  transfer  function.  Since  the  locations  of  the  sidebands 
are  determined  by  the  modulation  frequency  cos (2hDx/a j,) ,  the  minimum 
separation  between  the  apertures  must  be  D  =  2d  as  illustrated  in 
Figure  10,  in  order  for  the  signal  to  be  separated  completely  from 
the  zero  order  term.  However,  to  achieve  satisfactory  SNR,  we  should 
make  D  >  2d.  First  of  all,  the  dc  bias  term  for  a  low  contrast  input 
is  much  larger  than  the  modulated  signal.  Unless  the  image  spectrum 
is  translated  far  enough  away  from  the  dc  term  (by  increasing  0) , 
the  zero  order  side'lobes  would  contribute  significantly  to  the  fil¬ 
tered  image.  Moreover,  system  noise  such  as  grain  noise  drops  sharp¬ 
ly  as  we  move  away  from  the  zero  order.  Better  SNR  can  therefore  be 
obtained  with  the  use  of  higher  modulation  frequencies. 
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While  D  should  be  maximized  for  optimum  SNR,  d  should  be  mini¬ 
mized.  The  system  noise  (e.g.,  grain  noise)  within  the  filter  band¬ 
width  may  be  approximated  to  be  white.  The  amount  of  noise  power 
transmitted  through  the  filter  would  then  be  proportional  to  the  fil¬ 
ter  bandwidth.  Since  the  filter  bandwidth  is  designed  to  be  equiva¬ 
lent  to  the  image  bandwidth,  the  SNR  is  maximized  by  minimizing  d. 
However,  we  cannot  arbitrarily  increase  D  and  decrease  d.  With  the 
system  shown  in  Figure  9,  the  separation  D  is  limited  by  the  lens 
aperture.  The  maximum  separation  would  be  0  =  A  -  d  where  A  is  the 
diameter  of  the  lens.  Nevertheless,  we  shall  later  show  that  this 
limitation  can  be  eliminated  with  the  use  of  a  grating  interfero¬ 
meter.  Decreasing  d  on  the  other  hand  would  simultaneously  reduce 
the  resolution  and  efficiency  of  the  system.  With  most  applications, 
the  resolution  requirements  for  achieving  identification  are  not  very 
high  and  we  find  that  generally,  the  lost  in  efficiency  is  the  more 
detrimental  effect. 

The  use  of  very  small  apertures  may  require  an  overly  long  expo¬ 
sure  to  compensate  for  the  low  light  gathering  efficiency  of  the  sys¬ 
tem.  The  minimum  aperture  size  is  thus  determined  by  the  longest 
exposure  time  allowable  by  the  particular  application. 

4.2  GRATING  INTERFEROMETERS 

The  modulation  frequency  achievable  with  the  dual-aperture  system 
is  limited  by  the  numerical  aperture  of  the  lens.  To  overcome  this 
limitation,  an  imaging  system  utilizing  a  grating  interferometer  can 
be  used.  One  possible  arrangement  is  shown  in  Figure  11.  The  first 
grating  separates  the  input  light  field  into  two  parts  and  the  second 
grating  recombines  them.  The  gratings  are  volume  gratings  of  the 
same  spatial  frequency  constructed  with  dichromated  gelatin.  The 
first  grating  is  made  to  have  a  diffraction  efficiency  of  50^  while 
the  second  grating  approaches  100%.  Thus  aside  from  absorption, 
nearly  all  of  the  input  light  energy  is  directed  to  the  output  plane. 
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FIGURE  H.  IMAGING  A  SCATTERING  MEDIUM  WITH  A  GRATING  INTERFEROMETER 
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It  can  be  shown  that  at  plane  2=0,  stationary  fringes  with  spatial 
frequency  identical  to  that  of  the  gratings  are  formed  regardless  of 
the  phase  variation  of  the  input  light  field2-**.  Thus  stationary 
fringes  will  be  formed  by  both  the  unscattered  light  and  the  spatial¬ 
ly  incoherent  scattered  light.  In  order  to  allow  only  the  unscat¬ 
tered  light  to  form  stationary  fringes,  a  lateral  shear  of  the  image 
field  is  produced  by  moving  the  output  plane  to  z  =  iz  as  shown  in 
Figure  12.  The  output  intensity  distribution  can  be  approximated  to 
be 

2 

B2  +  A2  1  f ( x ,  y)  *  h(x,  y)| 

|[1  +  cos  (2* Gqx  ♦  e(x,  y)]j  (4-3) 

2  2 

where  B  and  A  are  the  intensities  of  the  scattered  and  unscat¬ 
tered  light  at  the  image  plane  respectively.  The  light  pattern  is 
thus  similar  to  that  obtained  with  the  dual-aperture  system  except 
that  the  modulation  frequency  is  now  determined  by  the  grating  fre¬ 
quency  Gq.  The  optimum  choice  of  Gq  is  dependent  on  the  charac¬ 
teristic  of  the  recording  material  utilized  with  the  system.  For 
many  materials  such  as  photographic  film,  the  amount  of  scattered 
noise  decreases  almost  exponentially  with  spatial  frequency.  On  the 
other  hand,  the  MTF  of  the  recording  material  also  decreases  with 
spatial  frequency.  The  fringe  contrast  would  be  reduced  if  the  modu¬ 
lation  frequency  is  too  high.  An  optimum  modulation  frequency  there¬ 
fore  exist  as  illustrated  in  Figure  13.  The  grating  frequency  Gq 
should  be  chosen  such  that  it  matches  the  characteristics  of  the  re¬ 
cording  material . 

To  process  the  recorded  transparency  to  extract  the  modulated 
target  image  from  the  background  bias,  a  coherent  optical  processor 
as  shown  in  Figure  14(a)  can  be  used.  The  Fourier  transformation  of 
the  input  transmittance  function  is  performed  by  the  converging 
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FIGURE  13.  DETERMINATION  OF  OPTIMUM  GRATING  FREQUENCY 


FIGURE  14B.  INCOHERENT  OPTICAL  PROCESSING  OF  TARGET  IMAGE 
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beam.  The  lens  aperture  is  placed  in  the  spatial  frequency  plane  at 
a  position  that  corresponds  to  the  carrier  frequency  of  the  modulated 
signal.  The  lens  aperture  would  act  as  the  bandpass  filter  and  the 
aperture  size  is  set  according  to  the  bandwidth  of  the  signal 
spectrum. 

The  input  transparency  is  imaged  onto  the  output  plane  by  the 
processing  lens.  Since  only  the  modulated  area  of  the  input  will 
diffract  light  to  the  lens  aperture,  the  modulated  target  image  will 
be  separated  from  the  zero  order  term  to  produce  an  output  image  with 
significantly  improved  contrast. 

The  recording  is  essentially  a  local  reference  beam  imaqe  plane 
hologram  and  as  with  other  types  of  image  plane  holograms,  incoherent 
light  can  be  used  to  reconstruct  the  holographic  image7.  One  ver¬ 
satile  approach  is  to  make  use  of  spatially  incoherent  light.  A 
light  source  with  varying  degrees  of  spatial  coherence  can  be  ob¬ 
tained  by  focusing  a  coherent  light  beam  on  a  rotating  ground  glass 
diffuser.  To  decrease  the  spatial  coherence,  one  may  simply  increase 
the  spot  size  of  the  focused  beam.  We  find  that  by  directly  illumi¬ 
nating  the  rotating  diffuser  with  a  laser  beam  as  shown  in  Figure 
14(b),  a  substantial  reduction  of  the  coherent  noise  can  be  achieved 
without  overly  affecting  the  SNR  of  the  processed  image.  Alterna¬ 
tively,  one  can  make  use  of  a  light  source  with  an  appropriate  spec¬ 
tral  bandwidth  to  create  similar  effects.  A  high  power  LED  would  be 
a  very  good  candidate  due  to  its  small  size  and  modest  cost. 

4.3  STABILITY  REQUIREMENTS 

The  major  advantage  offered  by  the  local  reference  beam  holo¬ 
graphic  technique  is  its  relatively  insensitivity  to  vibration  and 
atmospheric  turbulence.  First,  let  us  examine  the  stability  require¬ 
ments  for  recording  a  conventional  hologram. 
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In  Figure  15,  we  illustrate  the  basic  geometry  for  recording  a 
hologram  with  an  off-axis  reference  beam.  The  light  field  impinging 
on  the  recording  plane  can  be  written  as 

i  .  i f id (t  -  tt\  +  4R(x,y)l 

E(x,  y)  =  |ER(x,  y)| e  L  \  vl/  J 

I  i L (t  -  ^0(x,y)l 

*  |e0<x,  y)|e  L  \  V  0  J  (4-4) 

where  V-j  and  are  the  velocities  of  propagation  through  the 
two  optical  paths  (V  =  c/n  where  c  is  the  velocity  of  light  in  free 
space  and  n  is  the  refractive  index).  Recording  the  light  intensity 
on  a  photographic  plate,  the  resulting  transparency  would  have  a 
transmittance  of 


Ta(x,  y)  *  l.ER{x,  y)|2  ♦  |eq(x,  y)| 


+  2 |er (x,  y)|  | E  ( x,  y)| 


(4-5) 


'[(‘r<x-  »>  -  ‘o(x-  *»)*  {77  -  T^) 


In  order  to  record  the  fringe  pattern  with  the  highest  possible 
contrast,  the  fringe  pattern  should  remain  stationary  during  the  ex¬ 
posure.  Any  changes  in  the  phase  of  the  sinusoidal  interference 
pattern  will  have  a  smearing  effect  on  the  recording. 

Let  us  first  assume  that  the  recording  material  is  rigidly 
mounted  and  the  largest  permissible  fringe  movement  is  one-tenth  of 
a  period.  Any  changes  in  LR,  LQ,  Vj  or  Vj  will  produce  a 
shift  in  the  interference  pattern.  From  Eq.  4-5  we  can  see  that  to 
satisfy  the  stability  requirement,  it  would  be  necessary  that 
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If  the  refractive  indices  in  the  two  light  paths  are  both  equal  to 
1,  we  obtain 


1  ?  (LR  ■  Lo>  s  f  (4-7> 

Illuminating  with  a  laser  source  that  emits  at  0.6328  um,  we  have 

A ( Lp  -  Lq)  <  63.28  nm 

That  is,  the  change  in  the  path  difference  during  exposure  must  be 
less  than  63.28  nm. 

The  change  in  the  path  difference  can  be  caused  by  movement  (vi¬ 
bration)  of  the  target  or  the  beam  steering  optics.  To  avoid  a 
change  in  path  difference  of  more  than  63.28  nm,  the  components  must 
not  be  allowed  to  move  for  more  than  a  fraction  of  that  amount.  For 
example,  if  the  illuminating  beam  is  reflected  off  a  component  at  an 
angle  e  as  illustrated  in  Figure  16,  a  longitudinal  displacement  of 
d  for  the  component  will  produce  a  change  in  path  length  of 

aL  -  co?T  11  *  cos  29>  <4-8> 

For  e  =  15“,  the  maximum  permissible  displacement  for  the  component 
would  only  be  32.76  nm. 

We  now  assume  that  the  components  in  the  optical  arrangement  are 
stationary  during  the  exposure  but  the  refractive  indices  vary. 
Since  (LR  -  Lq)  <  coherent  length,  we  may  approximate  LR  = 
Lq  =  L  and  express  Eq.  (4-6)  as 

4L(vj  ‘  rj\ =  &[4(ni  -  n2}]  -  i  (4_9) 

For  L  -  1  km,  the  changes  in  the  refractive  indices  must  be  less  than 
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recording 


FIGURE  15.  OFF-AXIS  RECORDING  OF  CONVENTIONAL 

HOLOGRAM 


FIGURE  16.  PATH  LENGTH  VARIATION  DUE  TO  DISPLACEMENT  OF  AN  OPTICAL  COMPONLN 


A ( n -j  -  n2)  <  3.98  x  10~10 

Inside  a  controlled  laboratory  requirement  it  is  possible  to  meet 
these  stringent  requirements.  However,  under  actual  field  condi¬ 
tions,  it  is  unlikely  that  these  requirements  can  be  satisfied.  Com¬ 
ponents  are  mounted  independently  and  widely  apart.  Their  relative 
displacement  due  to  normal  vibration  can  easily  exceed  the  acceptable 
limit.  Moreover,  these  components  are  separated  by  a  dynamic  air 
mass.  The  refractive  indices  of  the  air  mass  changes  constantly  and 
with  a  variation  much  greater  than  3.98  x  10~^.  We  can  therefore 
conclude  that  the  conventional  hologram  would  not  be  practical  in 
actual  field  conditions. 

To  examine  the  stability  requirement  for  recording  a  local  refer¬ 
ence  beam  hologram  with  our  grating  interferometric  imaging  system, 
we  shall  first  assume  that  the  recording  material  is  rigidly  mounted 
and  free  of  vibrations.  Once  again,  to  ensure  that  the  fringes  would 
not  move  for  more  than  one-tenth  of  a  period,  it  is  necessary  that 

a(Lr  -  Lq)  S  63.28  nm 

With  the  grating  interferometer  however,  the  two  interfering 
light  beams  travel  through  a  nearly  identical  light  path.  For  exam¬ 
ple,  with  the  system  illustrated  in  Figure  17,  a  lateral  shear  of  ax 
between  the  two  interfering  images  would  correspond  to  a  shear  of  S 
=  1000  ax  at  the  target  plane  1  km  away.  If  the  target  is  under  the 
illumination  of  a  collimated  beam  incident  at  an  angle  e,  the  light 
paths  of  two  interfering  light  waves  may  be  illustrated  as  in  Figure 
18.  Displacing  the  target  longitudinally  by  a  distance  d  will  pro¬ 
duce  a  change  in  path  difference  of 
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For  example,  if  S  »  100  mm  (or  ax  =  0.1  mm)  and  R  =  1  km,  in  order 
to  satisfy  the  condition  that  a(Lr  -  L  )  <  63.28  nm,  the  longi¬ 
tudinal  motion  of  the  target  must  be  limited  to 

d  <  2.4  mm 


If  the  target  moves  laterally,  the  main  effect  is  to  produce  a  cor¬ 
responding  lateral  shift  of  the  fringes.  However,  due  to  the  demag¬ 
nification  factor  of  1000,  the  amount  of  permissible  lateral  motion 
would  be  d  <0.1  mm  for  a  1000  lines/mm  modulation  frequency. 

The  use  of  the  local  reference  beam  holographic  technique  with  a 
grating  interferometer  has  eliminated  the  need  of  separate  beam 
steering  optics  for  the  two  interfering  light  beams.  And  as  we  have 
shown,  it  also  relaxes  substantially  the  stability  requirement  for 
the  target  to  a  level  that  can  be  satisfied  in  the  field. 

We  had  assumed  that  the  recording  material  is  free  of  any  vibra¬ 
tion.  The  permissible  amount  of  motion  for  the  recording  material 
(and  the  receiving  optics)  is  the  same  for  the  conventional  and  local 
reference  beam  techniques,  that  is,  one-tenth  of  a  period  of  0.1  ym. 
Fortunately,  one  has  full  control  of  the  receiving  side  of  the  system 
and  motions  within  the  receiver  can  be  easily  reduced  to  the  accept¬ 
able  level  by  rigid  mounting  and  vibration  isolation. 

The  changes  in  the  refractive  index  of  the  air  mass  in  the  light 
path  have  a  severe  effect  on  the  conventional  holographic  technique 
because  the  interfering  light  beams  have  to  travel  through  very  dif¬ 
ferent  light  paths.  With  the  local  reference  beam  technique,  the 
interfering  beams  travel  through  a  nearly  identical  path.  The 
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changes  in  the  refractive  index  of  the  air  mass  would  affect  both 
beams  equally  and  thus  maintain  the  relative  phase  difference  between 
the  beams.  This  would  be  true  as  long  as  the  spatial  variation  of 
the  refractive  indices  in  the  air  mass  is  smaller  than  S,  the  separa¬ 
tion  of  the  interfering  point  radiators. 

There  are  also  other  advantages  in  utilizing  the  local  reference 
beam  technique.  In  order  for  two  light  beams  to  interfer,  their  path 
difference  must  be  less  than  the  coherent  length  of  the  illuminating 
light  source.  For  a  He-Ne  laser,  the  coherence  length  is  typically 
about  40  cm.  In  making  a  conventional  hologram,  a  precise  prior 
knowledge  of  the  target  range  would  therefore  be  necessary  such  that 
the  light  path  can  be  equalized.  Since  the  target  is  supposed  to  be 
obscured,  this  information  cannot  be  assumed  available.  Moreover, 
even  if  the  target  range  is  known,  the  depth  of  the  target  field  must 
be  restricted  to  that  of  the  coherent  length.  With  the  local  refer¬ 
ence  beam  technique,  the  path  lengths  of  the  two  interfering  beams 
would  automatically  equalize  since  they  are  generated  from  point  ra¬ 
diators  on  the  same  surface.  For  the  same  reason,  the  target  field 
can  be  of  any  depth.  For  most  applications,  a  light  source  with  a 
coherent  length  of  less  than  1  cm  would  provide  adequate  temporal 
coherence.  This  would  allow  a  much  greater  choice  of  light  sources, 
making  it  feasible  to  utilize  devices  such  as  laser  diodes  and  CO,, 
1 asers. 

4.4  PERFORMANCE  ESTIMATION 

To  determine  the  effectiveness  of  such  a  system  in  imaging 
through  fog,  we  examined  a  specific  system  as  illustrated  in  Figure 
19.  The  bandwidth  of  the  image  obtained  with  the  system  would  be 
2d/xL  =  250  lines/mm  and  the  modulation  frequency  is  equal  to  sin 
e/x  =  540  lines/mm.  The  optical  processing  arrangement  shown  earlier 
in  Figure  14(a)  is  used  to  perform  the  spatial  filtering.  The  zero 
order  and  the  modulation  terms  are  separated  by  the  diffraction  of 
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the  converging  beam.  Bandpass  filtering  is  accomplished  by  setting 
the  lens  aperture  according  to  the  Image  bandwidth.  To  assure  maxi¬ 
mum  fringe  contrast,  the  exposure  should  be  made  such  that  the  re¬ 
cording  is  biased  at  a  point  where  the  material  exhibits  the  highest 
transfer  characteristic  (i.e.,  the  region  with  the  steepest  slope  In 
the  T^-E  curve).  For  example,  with  Kodak  649F  plate  as  the  record¬ 
ing  material,  the  plate  should  be  exposed  to  give  an  average  trans¬ 
mittance  of  *  0.5.  The  percentage  of  input  light  diffracted  to 
one  of  the  sidebands  would  be  about 


TAe  H(w)A^"p 
2(A2  +  B2)_ 


x  100% 


(4-11) 


where  e  is  the  slope  of  the  transfer  characteristic  at  TA  and  H(w) 
is  the  frequency  response.  For  649  F  plates,  a  =  1.25  at  T^  «  0.5 
and  H(w)  >  0.8  at  u  *  500  lines/mm.  If  the  imaging  process  between 
the  input  transparency  and  the  output  plane  is  one-to-one,  then  the 
intensity  of  the  speckled  target  image  would  be  equal  to 


0.0625 


P  ? 
Tr  W/mm 


(4-12) 


where  P  is  the  input  light  power  and  S  is  the  total  area  of  the  input 
transparency  illuminated.  Although  this  is  only  a  first  order  ap¬ 
proximation,  it  provides  a  good  estimate  of  the  input  Image  intensity 

2  2 

when  B  »  A  (I.e.,  very  small  Input  modulation). 

The  main  source  of  noise  at  such  a  high  diffraction  angle  is  the 
scattering  noise  of  the  photographic  plate.  For  649F,  at  540 
lines/mm  and  T^  -  0.5,  the  power  of  the  scattering  noise  is  about 
3  x  10~9  P  W/line/mm2.  If  the  filter  bandwidth  is  made  to  be 
equal  to  the  signal  bandwidth,  that  is  250  lines/mm,  the  total  noise 
power  that  is  transmitted  to  the  output  plane  would  be 


1 


3  x  1(T9  x  « 


P  watts 


Since  the  noise  power  is  evenly  distributed  across  the  image  area  S, 
the  intensity  of  the  scattering  noise  would  be  equal  to 

I$  =  1.47  x  10"4  £ 

The  SNR  of  the  processed  image  would  therefore  be 
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corresponding  to  an  output  contrast  of 

0.0625 


‘'out 


2(1.47  x  10"4)  +  0.0625 


x  100  %  (4-14) 
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For  B  »  A  ,  we  may  make  the  following  approximation 


2  2  2 
A  _  A^  _  hc 
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The  output  contrast  CQUt  can  then  be  written  in  terms  of  input  con¬ 
trast  (as  observed  directly  through  the  scattering  medium)  as 


2.5  x  10-5  Cin2 


"out 


2.94  x  10"4  +  2.5  x  10-5  Cin2 


x  100% 


(4-15) 


where 


'in  =  77 


A  +  2B 


x  100% 


We  plotted  in  Figure  20  the  estimated  contrast  of  the  processed  image 
as  a  function  of  the  input  contrast.  We  find  that  the  processed 
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A  -  540  lines/mm  modulation,  filter 

bandwidth  equal  to  signal  bandwidth 


B  -  540  lines/mm  modulation,  filter 

bandwidth  equal  to  1/2  signal  bandwidth 


C  -  1000  lines/mm  modulation,  filter 

bandwidth  equal  to  1/2  signal  bandwidth 


- 1 - 1 — , — r- 

Input  Image  Contrast 


10 


FIGURE  20.  ESTIMATED  SYSTEM  PERFORMANCES 


image  would  be  recognizable  by  a  human  observer  if  the  input  image 
contrast  is  0.5%  or  higher.  Another  indicator  of  system  performance 
is  the  break  even  point.  A  positive  gain  in  image  contrast  would  be 
obtained  if  the  input  image  contrast  is  greater  than  0.12%. 

To  improve  the  performance  of  the  system,  the  simplest  way  is  to 
decrease  the  filter  bandwidth.  We  note  that  the  noise  power  is  even¬ 
ly  distributed  across  the  filter  passband  while  the  power  of  the 
image  spectrum  decreases  linearly  from  the  center  frequency.  Thus, 
if  we  decrease  the  filter  bandwidth  by  half,  the  noise  power  would 
be  cut  by  a  factor  of  4  while  the  signal  power  would  only  be  lowered 
by  a  factor  of  2.  The  resulting  SNR  would  effectively  be  doubled. 
The  performance  of  the  system  using  a  filter  bandwidth  half  the  size 
of  the  signal  spectrum  is  also  plotted  in  Figure  20.  The  processed 
image  will  be  recognizable  with  an  input  image  contrast  as  low  as 
0.35%  and  the  break  even  point  is  reduced  to  0.05%.  Reducing  the 
filter  bandwidth  would  also  produce  a  marginal  loss  in  image  resolu¬ 
tion.  Nevertheless,  the  available  resolution  would  still  be  higher 
than  needed  in  most  applications,  especially  if  the  output  is  to  be 
displayed  on  a  TV  screen.  The  marginal  loss  in  resolution  would 
generally  be  an  acceptable  trade-off  for  better  SNR. 

To  further  improve  the  SNR,  we  can  increase  the  modulation  fre¬ 
quency.  For  example,  in  Figure  20,  we  plotted  the  system  using  a 
1000  lines/mm  modulation  frequency.  A  human  observer  would  be  able 
to  recognize  the  processed  image  if  the  input  image  contrast  is  0.28% 
or  higher.  However,  we  cannot  continuously  obtain  a  SNR  improvement 
by  increasing  the  modulation  frequency.  The  finite  frequency  re¬ 
sponse  of  the  recording  material  would  place  a  limit  on  the  modula¬ 
tion  frequency  as  illustrated  earlier  in  Figure  13.  With  Kodak  649F 
as  the  recording  material,  the  optimum  modulation  frequency  is  around 
1000  lines/mm,  corresponding  to  a  diffraction  angle  of  40’  for  0.6328 
urn  illumination.  Other  factors  also  impose  additional  restraints  on 
the  modulation  frequency  of  the  grating.  The  grating  became  more 
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angularly  selective  with  frequency.  That  is,  the  diffraction  effi¬ 
ciency  falls  quickly  as  the  incoming  beam  moves  off  the  Bragg  angle. 
This  would  produce  a  decrease  in  the  intensity  of  image  points  away 
from  the  optical  center,  resulting  in  a  reduction  of  the  field  of 
view.  Furthermore,  at  high  diffraction  angles,  the  off-axis  aberra¬ 
tions  are  more  severe.  The  problem  of  aberrations  will  be  discussed 
in  a  later  section. 

The  use  of  the  local  reference  beam  technique  makes  the  holo¬ 
graphic  system  practical  in  field  conditions.  However,  a  price  in 
SNR  is  paid  in  gaining  these  advantages.  Comparing  Eq.  3-3  and  Eq. 
4-3,  we  see  that  with  the  conventional  holographic  technique,  |f(x, 
y)  *  h(x,  y)|  is  reconstructed.  While,  with  the  local  reference  beam 
technique,  [f(x,  y)  *  h(x,  y)| c  is  reconstructed.  Since  the  band- 

r\ 

width  of  ff(x,  y)  *  h(x,  y)|  L  is  twice  as  wide  as|f(x,  y)  *  h(x, 
y)|,  the  SNR  would  be  a  factor  of  4  lower  for  the  local  reference  beam 
technique.  If  the  filter  bandwidth  is  set  to  be  half  of  the  signal 
bandwidth,  then  the  SNR  loss  would  be  reduced  to  a  factor  of  2. 
Thus,  the  conventional  holographic  technique  should  be  used  whenever 
the  experimental  conditions  allow 

4.5  RECORDING  MATERIAL 

The  quantitative  analysis  presented  in  the  previous  section  ap¬ 
plies  only  to  the  chosen  recording  material,  the  Kodak  649F  spectro¬ 
scopic  plate.  The  649F  plate  was  chosen  because  of  its  common  usage 
in  holography.  Its  performance  would  serve  well  as  a  standard  with 
which  other  materials  may  be  compared. 

Upon  choosing  a  recording  material  for  the  interferometric  imag¬ 
ing  system,  the  optimum  modulation  frequency  and  optimum  exposure 
must  then  be  determined.  The  optimum  modulation  frequency  is  ob¬ 
tained  by  measuring  the  SNR  as  a  function  of  spatial  frequency.  For 
649F,  the  optimum  modulation  frequency  is  around  1000  lines/mm  as 
shown  in  Figure  21.  The  optimum  exposure  can  be  similarly  determined 
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FIGURE  21.  CHARACTERISTICS  OF  649F  PLATES  AS  A  FUNCTION  OF  SPATIAL  FREQUENCY 


by  measuring  the  SNR  as  a  function  of  exposure.  The  behavior  of  649F 

as  a  function  of  exposure  under  0.6328  ym  illumination  is  plotted  in 

2 

Figure  22.  The  optimum  exposure  Is  achieved  at  about  100  yj/cm  . 
It  should  be  noted  that  the  SNR  measurements  are  dependent  on  the 
modulation  depth  of  the  input  signal.  For  an  input  with  high  modu¬ 
lation,  the  film  nonlinearity  will  produce  intermodulation  noise  that 
cause  the  optimum  exposure  to  occur  at  a  lower  exposure  level.  In 
order  to  obtain  optimum  parameters  appropriate  for  our  present  appli¬ 
cation,  the  SNR  measurements  were  made  using  a  fringe  pattern  with 
fairly  low  modulation  (282). 

One  of  the  most  serious  drawbacks  of  using  649F  is  its  low  sen¬ 
sitivity.  In  situations  where  the  input  light  intensity  is  very  low, 
the  required  exposure  time  would  be  impractical ly  long.  Various 
techniques  may  be  used  to  increase  the  speed  of  the  emulsion.  For 
example,  by  soaking  the  photographic  plate  in  water  for  about  15 
minutes  inside  a  liquid  gate  before  the  exposure,  the  sensitivity  of 
the  emulsion  can  be  increased  by  as  much  as  4  times.  The  use  of 
special  development  procedure  may  increase  the  speed  by  another  fac¬ 
tor  of  2.  However,  to  improve  the  system  sensitivity  significantly, 
the  substitution  with  a  faster  recording  material  would  be  necessary. 

Another  fine  grain  holographic  recording  material  is  the  Kodak 

131  plate.  The  optimum  exposure  for  the  131  plate  is  about  500 

2 

nJ/cm  ,  making  it  200  times  faster  than  the  649F  plate.  Another 

possible  substitution  is  the  Kodak  AHU5460  microfilm.  The  film  has 

a  low  noise  acetate  substrate  and  it  has  been  utilized  for  many  years 

in  the  recording  of  data  for  optical  processing  of  synthetic  aperture 

radar  (SAR)  signals.  The  optimum  exposure  for  the  5460  microfilm  is 
2 

about  100  nJ/cm  .  However,  to  make  a  photographic  emulsion  more 
sensitivity,  it  is  necessary  to  use  larger  silver  halide  crystals. 
The  larger  grain  size  would  produce  more  severe  grain  noise.  Thus, 
a  tradeoff  between  sensitivity  and  SNR  exists  and  the  choice  of  the 
recording  material  would  depend  on  the  requirements  of  the  specific 
application. 
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Kodak  649F  developed  in  D-19 


Like  all  silver  halide  photographic  material,  649F  requires 
chemical  developments.  The  normal  development  process  requires  about 
8  minutes  (5  minutes  development  and  3  minutes  fixation).  The  de¬ 
velopment  time  can  be  shortened  by  the  use  of  a  higher  development 
temperature,  a  faster  acting  developer  or  a  higher  concentration  of 
hydroxide  in  the  developer.  It  is  possible  to  reduce  the  processing 
time  to  less  than  a  minute.  However,  rapid  development  usually  re¬ 
sults  in  a  lower  quality  recording.  The  developed  silver  grains  tend 
to  clump  together  and  concentrate  near  the  emulsion  surface.  The 
result  is  a  higher  scattering  noise  and  poorer  uniformity. 

To  substantially  reduce  the  processing  time,  it  would  be  neces¬ 
sary  to  use  a  real  time  light  modulator  in  place  of  the  photographic 

material.  Unfortunately,  most  real  time  light  modulators  such  as 

9  10 

PROM  and  liquid  crystal  light  valves  are  sensitive  predomi¬ 
nantly  in  the  blue-UV  region.  The  photoplastic  device  has  been  shown 
to  be  sensitive  to  near  IR  illumination,  responding  up  to  1.15  vm 
and  its  sensitivity  is  compatible  with  that  of  the  649F  plate14. 
However,  the  achievable  SNR  is  lower15. 

When  a  photographic  material  is  used  with  the  system,  the  optical 
processing  should  be  performed  with  the  plate  or  film  inside  a  liquid 
gate.  The  emulsion  surface  is  not  flat  and  the  random  thickness 
variation  across  the  emulsion  would  introduce  additional  scattering 
noise.  By  placing  the  film  inside  a  liquid  whose  refractive  index 
matches  that  of  the  emulsion,  the  effect  of  the  irregular  surface 
relief  can  be  eliminated.  Xylene  with  a  refractive  index  of  1.5  is 
the  most  commonly  used  liquid.  It  matches  particularly  well  acetate 
based  film  such  as  the  5460  microfilm. 

For  the  best  combination  of  sensitivity,  processing  speed  and 
SNR  the  photographic  plate  or  film  can  be  placed  inside  a  liquid  gate 
filled  with  water.  The  film  is  exposed  while  it  is  inside  the  water 
gate  and  the  film  is  then  taken  out  for  development.  (It  is  also 
possible  to  develop  the  film  in  situ  by  replacing  the  water  with 
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developer.)  After  a  brief  fixing  period  (just  enough  to  clear  the 
film),  the  developed  transparency  is  placed  right  back  into  the  water 
gate  for  optical  processing.  Placing  the  film  inside  the  water  gate 
before  the  exposure  increases  the  film  sensitivity.  Presoaking  the 
emulsion  would  also  facilitate  the  penetration  of  the  developer,  re¬ 
sulting  in  a  faster  and  more  complete  development.  Putting  the  de¬ 
veloped  film  right  back  into  the  water  gate  eliminates  the  drying 
step.  Although  water,  with  a  refractive  index  of  1.3,  does  not  match 
the  emulsion  as  well  as  xylene,  a  substantial  amount  of  the  surface 
relief  effect  can  still  be  eliminated. 

4.6  IMAGE  ABERRATIONS  AND  DESIGN  PARAMETERS 

In  this  section,  we  shall  examine  the  design  parameters  for  a 
grating  interferometer  that  will  produce  minimum  image  aberrations. 
First,  let  us  consider  the  case  where  a  sinusoidal  plane  grating  is 
placed  normal  to  the  optical  axis  of  the  imaging  system  and  examine 
how  well  it  can  image  an  object  point  located  at  the  optical  center. 
One  approach  is  to  use  the  ray  tracing  technique  to  determine  the 
shape  and  location  of  the  image  formed  by  the  diffracted  light. 

A  light  ray  with  directional  cosines  cos  ex  and  cos  ey  impinging 
on  the  grating  will  emerge  as  3  rays,  the  ±  first  order  and  the  zero 
order.  Only  one  of  the  first  order  diffracted  rays  is  of  interest 
to  us.  The  directional  cosines  of  this  diffracted  ray  can  be  deter¬ 
mined  by 

cos  ex  =  cos  tx  -  ^ 
and 


cos  ey 


cos  e 

y 
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where  dx  and  dy  are  the  fringe  spacings  in  the  corresponding 
direction. 

In  Figure  23(a),  we  show  the  ray  trace  for  the  point  image  formed 
by  the  diffracted  light.  We  see  that  the  point  image  is  severely 
aberrated.  In  Figure  23(b),  we  show  the  experimentally  obtained 
light  distributions  at  three  different  image  planes.  Coma  and  astig¬ 
matic  aberration  is  quite  evident.  An  alternate  approach  in  examin¬ 
ing  the  image  aberration  is  to  consider  the  grating  as  a  hologram 
constructed  with  two  point  sources  located  at  infinity  (i.e.,  plane 
waves).  The  hologram  is  then  reconstructed  by  a  converging  beam  that 
focuses  at  a  distance  Rc  from  the  plate.  If  we  confine  the  refer¬ 
ence  beam,  the  object  beam,  the  reconstructing  beam  and  the  image 
beam  all  in  the  x-z  plane,  the  reconstructed  point  image  is  found  to 
be  located  at  a  distance  fT  from  the  center  of  the  plate  as  illus¬ 
trated  in  Figure  24,  where  is  determined  by 

rrk'ik'ti 

Since  RQ  =  Rr  =  °°  ,  we  have  R^  =  Rc. 

The  wavefront  derivation  (in  wavelengths)  is  computed  by 
Champagne^  ^  to  be 

A  (x  +  r  >s 

+  \  (x2  +  y2)*cx 

- I (x2  v] 

where 


(spherical  aberration) 


(4-19) 
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sin  a  /sin  a  sin  a  \  Sin  a . 

Cx  =  -z-1  -  [— y-a  -  —?A - r2  (4-20) 


c 

\coma) 


.2  /  .  2  •  2  \  .2 
sin  a  /sin  a.  sin  a\  sin  a, 

Ax  =  -y— S-I-tt-S - - *-1  (4-21) 


c  \  Ro  Rr 

(astigmatic  aberration) 


Ri 


with  Rn  =  R 
o  r 


=  °°,  they  can  be  simplified  to 

c  1  1 

"I  R? 


(4-22) 


sin  a  sin  a. 

c  -  c _ I 

Sc  -  R. 


(4-23) 


and 


Ax  = 


.2  .2 

Sin  ac  Sin 

~c  "T 


(4-24) 


Since  =  Rc,  we  have  S  =  0;  there  would  not  be  any  spherical 
aberration.  By  placing  the  grating  perpendicular  to  the  optical  axis 
of  the  reconstructing  beam,  ac  «  0,  and  we  have 

sin  a.  sin2  a. 

C  =  - 1  and  A  =  - 1  (4-25) 

*  R.  R. 

The  image  would  therefore  suffer  from  coma  and  astigmatic 
aberration.  We  can  see  that  the  aberration  increases  with  the  dif¬ 
fraction  angle  and  astigmatism  is  the  more  dominant  aberration. 

From  Eqs.  4-23  and  4-24,  we  also  find  that  the  aberrations  can 
be  made  to  vanish  by  making  ac  =  a^.  This  can  be  achieved  by 
tilting  the  grating  an  angle  a./2.  In  Figure  25(a),  we  show  the 
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ray  trace  with  the  tilted  grating  arrangement.  The  aberrations  in¬ 
deed  disappear. 

What  we  have  described  would  apply  directly  to  the  lower  branch 
of  the  grating  interferometer.  The  upper  branch  on  the  other  hand, 
is  diffracted  twice.  By  tilting  both  gratings  by  an  angle  <V2, 
the  aberrations  will  be  corrected  as  Illustrated  In  Figure  25(b). 
Since  the  required  angle  of  tilt  for  the  second  grating  is  the  same 
with  the  upper  and  lower  branches  of  the  interferometer.  A  single 
solid  grating  can  therefore  be  used.  In  order  for  the  Images  re¬ 
constructed  by  both  branches  to  coincide,  the  central  rays  should 
form  a  parallelogram.  Together  with  the  angle  requirements,  an 
aberration-free  grating  Interferometer  can  be  designed  with  two  sinu¬ 
soidal  gratings  as  shown  in  Figure  26. 

The  use  of  plane  gratings  is  not  very  efficient  since  one-third 
of  the  light  is  diffracted  to  the  conjugate  order.  The  use  of  volume 
gratings  would  confine  the  diffraction  to  a  single  order.  By  making 
the  first  grating  to  have  a  diffraction  efficiency  of  50%,  and  the 
second  to  have  a  diffraction  efficiency  of  100%,  all  the  transmitted 
light  energy  would  be  directed  to  the  image  plane.  The  tilted  grat¬ 
ing  arrangement  can  also  be  used  with  volume  gratings  to  produce 
aberration-free  Images.  However,  if  the  grating  is  constructed  with 
two  plane  waves,  the  bragg  condition  would  be  satisfied  fully  by  only 
the  central  rays  as  illustrated  in  Figure  27(a).  The  peripheral  rays 
would  be  diffracted  less  efficiently.  The  result  is  a  reduction  in 
the  effective  aperture  of  the  imaging  system  in  the  direction  of  dif¬ 
fraction.  This  effect  is  particularly  pronounced  when  a  high  diffrac¬ 
tion  angle  is  used  with  a  fast  lens.  To  avoid  this  effect,  a  volume 
grating  that  is  constructed  with  two  point  sources  instead  two  plane 
waves  should  be  used.  The  volume  grating  is  placed  with  its  back  to 
the  impinging  converging  light  beam  as  Illustrated  In  Figure  27(b). 
The  bragg  condition  will  be  satisfied  across  the  entire  aperture. 
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FIGURE  27A.  DIFFRACTION  OF  A  CONVERGING  BEAM 
WITH  A  GRATING  CONSTRUCTED  WITH  PLANE  WAVES. 


FIGURE  27B.  DIFFRACTION  OF  A  CONVERGING  BEAM  WITH 
A  GRATING  CONSTRUCTED  WITH  TWO  DIVERGING  BEAMS. 


With  these  preliminary  discussions,  we  can  now  proceed  to  de¬ 
scribe  the  optical  arrangement  that  was  used  to  fabricate  the  grat¬ 
ings  utilized  in  our  interferometer. 

4.7  GRATING  FABRICATIONS 

The  arrangement  for  fabricating  the  first  grating  is  rather  sim¬ 
ple.  Two  point  sources,  placed  symmetrically  with  respect  to  the 
holographic  plate  are  located  at  a  distance  L  from  the  center  of  the 
plate.  L  is  determined  by  the  desired  grating  position  in  the  inter¬ 
ferometer  as  shown  in  Figure  28  (a).  When  the  developed  holographic 
grating  is  placed  at  the  appropriate  position  behind  the  imaging 
lens,  the  converging  beam  will  be  diffracted  to  form  an  unaberrated 
point  image  as  shown  in  Figure  28(b). 

To  generate  the  second  grating,  we  make  use  of  the  same  two  point 
sources  as  the  first  grating,  and  add  a  third  point  source  mid-way 
between  them.  The  location  of  this  third  point  source  will  determine 
the  eventual  image  position  in  the  grating  interferometer.  The  holo¬ 
graphic  plate  is  placed  midway  between  the  first  grating  and  the 
image  plane  such  that  the  central  rays  of  the  interferometer  will 
form  a  parallelogram.  The  light  emanated  from  this  third  point 
source  must  have  a  wide  enough  divergence  to  cover  the  entire  holo¬ 
graphic  plate.  This  would  require  a  lens  with  an  impossibly  large 
numerical  aperture.  To  overcome  this  problem,  we  use  a  lens  with  a 
numerical  aperture  slightly  larger  than  that  of  the  imaging  lens  and 
focus  a  beam  onto  a  volume  grating.  This  volume  grating  should  have 
a  spatial  frequency  identical  to  that  of  the  first  grating  and  a  dif¬ 
fraction  efficiency  of  50  %  (or  we  may  simply  make  use  of  the  first 
grating  that  was  generated).  The  undiffracted  beam  will  diverge  and 
illuminate  the  region  of  the  plate  where  the  upper  point  source  is 
also  illuminating.  The  diffracted  light  on  the  other  hand,  will  be 
directed  towards  the  lower  portion  of  the  plate  that  is  illuminated 
by  the  lower  point  source  as  illustrated  in  Figure  30(a).  Recorded 


FIGURE  28A.  FABRICATION  OF  FIRST  GRATING 
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FIGURE  28B.  DIFFRACTION  OF  A  CONVERGING  BEAM  BY  FIRST  GRATING 


FIG'.r:  29  A.  FAAPIC/'TION  OF  SECONO  GRATING 


FIGURE  29B.  IMAGING  GRATING  INTERFEROMETER 


on  the  holographic  plate  will  be  two  separate  sets  of  fringe  pat¬ 
terns.  The  spatial  frequencies  of  these  two  patterns  are  the  same 
but  the  fringes  are  tilted  differently.  After  the  plate  de¬ 
veloped,  we  can  assemble  with  the  two  gratings  an  interferometer  as 
shown  in  Figure  29(b).  The  interferometric  system  will  be  able  to 
produce  the  image  of  an  on-axis  point  object  without  any  aberration. 

To  generate  the  gratings,  dichromated  gelatin  was  used  as  the 

recording  material.  Dichromated  gelatin  is  the  best  holographic  re- 
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cording  material  available  today  .  The  maximum  achievable  dif¬ 
fraction  efficiency  is  over  90%  with  most  of  the  optical  loss  due  to 
absorption.  Since  the  gelatin  is  grainless,  the  material  also  has 
the  advantage  of  low  noise.  The  scattering  that  exists  witn  dichro¬ 
mated  gelatin  is  due  to  microscopic  cracks  in  the  gelatin  caused  by 
the  dehydration  operation.  Nevertheless,  the  noise  behavior  is  supe¬ 
rior  even  to  an  unbleached  649F  plate. 

The  gelatin  layer  can  be  coated  on  a  glass  substrate  by  dip¬ 
coating.  However,  photographic  materials  with  a  thick  emulsion  offer 
a  convenient  source  of  high  quality  gelatin  coatings.  The  familiar 
Kodak  649F  plate  was  used  to  generate  our  gratings.  The  gelatin 
layer  is  prepared  by  first  dissolving  all  the  silver  halide  crystals 
in  the  gelatin  with  a  fixer.  The  plate  is  then  conditioned  and  di¬ 
chromated.  After  the  exposure,  the  dichromated  gelatin  plate  is  de¬ 
veloped  by  dehydration  with  alcohol  baths.  The  preparation  and  de¬ 
velopment  procedures  for  dichromated  gelatin  plates  are  listed  in 
Tables  4-1  and  4-2. 

The  diffraction  efficiency  of  the  gratings  can  be  controlled  by 
varying  the  exposure  or  K-ratio  in  the  recording.  The  first  grating 
in  our  interferometer  both  diffracts  and  transmits  40  %  of  the  input 
light.  The  rest  is  lost  due  to  absorption,  scattering,  and  reflec¬ 
tion.  The  second  grating  has  diffraction  efficiency  of  75%.  Thus, 
the  overall  efficiency  of  the  grating  interferometer  is  about  60%. 


TABLE  4-1 

PREPARATION  OF  DICHROMATED  GELATIN  PLATES 
BASED  ON  649F 

DIRECTION 

1.  Soak  in  fixer  without  hardener  for  10  minutes  Room  light 

2.  Wash  in  running  water  at  90 “F  for  5  minutes;  start 
at  68“F  and  raise  temperature  gradually  to  90°F; 
following  wash,  lower  temperature  at  same  rate  to 

68° F.  Room  light 

3.  Soak  in  fixer  with  3.25%  hardener  for  10  minutes  Room  light 

4.  Wash  in  running  water  for  10  minutes  Room  light 

5.  Rinse  in  distilled  water  for  5  minutes  Room  light 

6.  Dip  into  a  Photo-Flo  200  solution  (2  drop/liter)  for 

30  seconds  and  remove  excess  with  sponge  Room  light 

7.  With  plate  in  horizontal  position,  dry  overnight 

at  room  temperature  Room  light 

8.  Soak  in  5%  ammoni urn- di chromate  solution  with 

2  drops/liter  of  photo-flow  200  for  5  minutes  Red  light 

9.  With  plate  in  horizontal  position,  dry  overnight 

at  room  temperature  Red  light 
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TABLE  4-2 

STANDARD  DEVELOPMENT  OF  DICHROMATED  GELATIN  PLATES 


DIRECTION 


1. 

Soak  in  0.5%  solution  of  ammonium  dichromate  for 

5  minutes 

Red 

light 

2. 

Soak  in  fixer  with  3. 25  %  hardener  for  5  minutes 

Red 

light 

3. 

Wash  in  running  water  for  10  minutes 

Room 

light 

4. 

Soak  in  distilled  water  for  2  minutes  with 
agitation 

Room 

1  i  ght 

5. 

Dehydrate  in  a  50/50  solution  of  distilled 
water  and  isopropanol  for  3  minutes 

Room 

1  i  ght 

6. 

Dehydrate  in  100%  isopropanol  for  3  minutes 

Room 

light 

7. 

Free-air  dry  for  at  least  1  hour 

Room 

1  i  ght 

8. 

Remove  water  completely  from  emulsion  by  baking 
or  with  a  vacuum  chamber 

Room 

1  i  ght 
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4.8  WIDEF I EL D  GRATING  INTERFEROMETER 

A  system  fabricated  with  the  technique  we  just  described  would 
be  capable  of  imaging  any  object  point  on  the  y-axis  perfectly.  How¬ 
ever,  for  the  object  points  that  are  located  away  from  the  y-axis, 
the  image  would  be  aberrated. 

We  have  shown  that  all  the  aberrations  can  be  eliminated  If  we 
make  a •  =  ac  and  we  designed  a  system  that  will  provide  aberra¬ 
tion  free  performance  for  an  object  point  at  the  optical  center. 
The  image  will  remain  unaberrated  as  the  object  point  moves  away  from 
the  optical  center  along  the  y-axis.  The  condition  that  -  <*c 
will  continue  to  be  satisfied.  If  the  object  point  moves  off  the 
y-axis  however,  the  point  image  will  suffer  from  coma  and  astigmatic 
aberration  as  illustrated  in  Figure  30  .  The  amount  of  aberrations 
would  be  proportional  to 


and 


Sin  (ac  +  Aa) 


sin  (a.j-  Aa) 


For  coma 


(4-26) 


.2  .2 
Sin  (ac  +  Aa)  Sin  (a,.  -  Aa) 


For  astigmatic  aberration 

where  Aa  is  the  corresponding  entrance  angle  of  the  off-set  object 
point  in  the  x-direction.  We  can  see  that  the  aberrations  increase 
with  aa.  Thus,  the  image  points  at  the  edge  of  the  image  field  would 
be  more  severely  aberrated  than  those  near  the  center  of  the  field. 
Furthermore,  we  find  that  coma  decreases  with  ac  while  astigmatism 
increases  with  ac  for  a  given  Aa.  The  effect  of  astigmatic  aberra¬ 
tion  nevertheless  dominates  and  a  grating  with  a  higher  spatial  fre¬ 
quency  would  be  more  sensitive  to  Aa. 
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FIGURE  30.  ABERRATION  OF  AN  OFF  AXIS  DIFFRACTION 


In  most  Army  applications,  the  grating  interferometer  would  be 
utilized  with  a  telephoto  Imaging  system.  The  field  of  view  (or 
field  angle)  would  generally  be  limited  and  the  off-axis  aberration 
would  not  be  significant  enough  to  be  a  problem.  However,  if  the 
application  demands  a  large  field  angle,  a  wide  field  grating  Inter¬ 
ferometer  can  be  constructed  with  the  use  of  computer  holographic 
techniques.  In  the  following,  we  shall  briefly  describe  the  concept 
and  fabrication  procedure  for  a  wide  field  grating  interferometer. 

For  each  field  angle,  we  can  compute  the  necessary  phase  trans¬ 
formation  t>n{p,  e)  to  produce  an  unaberrated  image.  For  example, 
for  an  object  point  «n,  we  can  propagate  the  light  wave  to  the 
plane  where  the  hologram  is  to  be  located  and  compute  phase  distri¬ 
bution  of  the  wavefront  in(p,  e)  at  that  plane.  Similarly,  we 
can  propagate  the  light  wave  (backward)  from  the  corresponding  image 
point  to  the  hologram  plane  and  obtain  the  necessary  wavefront 
<>n  out(p,  t0  ProdlJce  an  unaberrated  image  point.  The  required 
phase  transformation  for  producing  a  point  image  at  «n  would  be 

*n(p’  e)  *  *n,out{p’  e)  ~  0n,in{p’  •>  f4-28* 

After  computing  the  phase  transformation  0p(p,  e)  Tor  a  represen¬ 
tative  set  of  object  points  in  the  field  of  view,  we  can  then  obtain 
an  average  phase  transformation  function 

N  . 

*(p.  [  $n(p.  «)  (4-29) 

n*l 

It  has  been  shown  that  a  hologram  that  provides  such  a  phase  trans¬ 
form  will  produce  the  smallest  mean  squared  phase  error  over  the  en¬ 
tire  field  of  view  [20],  The  image  point  reconstructed  would  be 
slightly  larger  than  a  diffraction  limited  spot  but  the  small  spot 
size  can  be  maintained  over  a  much  larger  range  of  field  angles^7. 


To  produce  the  necessary  wavefront  to  construct  a  hologram  that 
gives  the  desired  phase  transformation  0(a,  e),  a  computer  generated 
hologram  can  be  used^®.  A  computer  generated  hologram  (CGH)  re¬ 
presents  the  amplitude  and  phase  of  a  wavefront  at  a  particular  sam¬ 
pled  position  by  the  size  and  position  of  a  binary  cell  within  a 
grid.  However,  the  reconstruction  efficiency  of  a  CGH  is  quite  low. 
And  due  to  the  limited  space-bandwidth  available  with  a  computer  dis¬ 
play,  the  carrier  frequency  and  numerical  aperture  of  the  CGH  have 
to  be  severely  restricted.  Thus,  instead  of  using  the  CGH  itself  in 
the  grating  interferometer,  it  is  copied  onto  a  dichromated  gelatin 
plate  as  illustrated  in  Figure  32.  We  note  that  with  such  a  scheme, 
the  CGH  is  not  constructed  to  produce  the  necessary  wavefront  0(p, 
9)  directly  at  the  hologram  plane,  but  it  is  made  to  produce  a  wave- 
front  which  upon  propagating  through  the  optical  system,  would  become 
«fi(p,  e)  at  the  recording  plane^. 

4.9  EXPERIMENTAL  RESULTS 

Before  the  system  was  evaluated  for  imaging  through  a  scattering 
medium,  it  was  tested  for  its  imaging  ability  using  a  resolution  tar¬ 
get  as  the  object.  The  target  transparency  was  diffusely  illuminated 
from  the  back  with  a  ground  glass  diffuser.  In  Figure  33(a),  we  show 
the  input  used  in  the  evaluation  and  in  Figure  33(b),  we  have  the 
image  obtained  with  the  grating  interferometric  system.  The  image 
quality  is  very  good  with  only  a  slight  degradation  at  the  edge  of 
the  field. 

The  image  in  Figure  33(b)  was  obtained  at  the  plane  where  the 
two  interfering  images  exactly  coincide  (i.e.,  z  =  0  plane).  In 
order  to  prevent  the  scattered  light  from  forming  stationary  fringes 
the  two  images  have  to  be  slightly  displaced  or  sheared.  To  obtain 
the  correct  amount  of  shearing,  a  back  illuminated  rotating  diffuser 
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FIGURE  31.  COMPUTATION  OF  REQUIRED  PHASE  TRANSFORMATION 

<P  (6ae) 


FIGURE  32.  FABRICATION  OF  WIDE-FIELD  DIFFRACTION  GRATING 
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was  inserted  between  the  target  plane  and  the  imaging  lens.  When 
the  gratings  are  perfectly  aligned,  stationary  fringes  can  be  ob¬ 
served  across  the  entire  image  field  with  a  microscope.  While  ob¬ 
serving  the  output  light  field,  the  second  grating  was  displaced 
sligntly  in  the  z-direction  until  the  fringes  disappeared.  Any  fur¬ 
ther  displacement  of  the  grating  would  result  in  an  excessive  image 
blur. 

To  evaluate  the  system  capability  in  imaging  through  scattering 
media,  a  qualitative  test  was  first  performed.  A  diffusely  illumi¬ 
nated  photographic  transparency  of  a  toy  tank  on  a  grass  lawn  was 
used  as  the  target.  To  produce  an  adjustable  scattering  medium,  a 
liquid  crystal  scattering  plate  was  used.  The  scattering  plate  was 
constructed  with  a  liquid  crystal  layer  sandwiched  between  two  glass 
plates.  The  liquid  crystal  was  made  to  operate  in  the  dynamic  scat¬ 
tering  mode  with  the  application  of  an  AC  voltage.  The  amount  of 
scattering  can  be  made  to  vary  by  adjusting  the  excitation  voltage. 
In  order  to  create  a  rapid  change  in  the  scattering  field,  the  scat¬ 
tering  plate  was  mounted  on  a  translator  stage  and  translated  across 
the  light  path  during  the  exposure. 

The  excitation  voltage  was  initially  set  at  zero  such  that  the 
liquid  crystal  plate  was  completely  clear  and  image  distinctly  vis¬ 
ible.  The  excitation  voltage  was  then  steadily  increased  until  the 
target  image  contrast  was  so  low  that  it  was  no  longer  recognizable. 
The  target  image  was  then  recorded  with  the  grating  interferometric 
system  and  processed  as  described  earlier.  In  Figure  34(a),  we  show 
the  target  image  as  observed  directly  through  the  scattering  medium 
and  in  Figure  34(b)  we  show  the  processed  image  obtained  with  the 

interferometric  system.  We  can  see  that  the  image  contrast  is  signi- 

» 

ficantly  improved  and  the  processed  target  image  is  easily 
recognizable . 

To  test  the  system  under  a  more  realistic  simulated  scattering 
environment,  a  glass  water  tank  filled  with  diluted  milk  was  used  as 
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FIGURE  34.  QUALITATIVE  T !  T  \,r\\  COMPLEX  TARGET  USING 
LIQUID  CRYSTAL  SCATTERING  P; AT!  AS  SCATTERING  MEDIUM 


the  scattering  medium.  The  suspension  of  the  fat  particles  in  a  3- 
dimensional  volume  creates  an  environment  very  similar  to  that  of 
natural  fog.  For  this  experiment,  a  toy  tank,  placed  just  outside 
the  glass  water  tank  and  illuminated  from  the  side,  was  used  as  the 
target.  The  concentration  of  the  milk  was  gradually  increased  until 
the  target  was  no  longer  recognizable.  The  target  image  was  then 
recorded  with  the  grating  interferometric  system.  The  processed 
image  is  shown  in  Figure  35.  The  target  is  easily  recognizable  as 
an  army  tank. 

It  is  difficult  to  define  the  image  contrast  of  a  complex  con¬ 
tinuous  tone  target.  Thus,  to  facilitate  a  quantitative  evaluation 
of  the  system,  we  used  a  simple  two-tone  object  as  illustrated  in 
Figure  36.  When  the  input  image  contrast  is  very  low,  I  and 
Imin  would  be  very  close  in  magnitude  and  it  would  be  difficult  to 
measure  the  input  contrast  accurately.  To  circumvent  this  problem, 
we  simulated  the  effect  of  imaging  through  a  scattering  medium  with 
the  arrangement  depicted  in  Figure  36.  With  such  an  arrangement, 
Imin  (scattered  light  intensity)  and  ^ax^min  (unscattered 
light  intensity)  can  be  measured  separately  by  blocking  the  other 
beam.  The  relative  intensities  of  signal  and  background  light 
fields  can  be  set  to  any  ratio  desired  by  adjusting  the  variable 
beam  splitters. 

To  evaluate  the  system  performance,  the  contrast  of  the  input 
image  was  adjusted  to  5%,  1%,  0.6%,  and  0.4%.  The  image  with  5$  con¬ 
trast  was  barely  recognizable  and  the  Images  with  the  lower  contrasts 
were  not  visible  at  all.  The  processed  image  obtained  with  the 
interferometric  system  were  8!$,  19%, 8.5%,  and  3%  respectively.  The 
experimental  results  are  plotted  in  Figure  37  together  with  the  esti¬ 
mated  performance  computed  earlier.  To  illustrate  the  clarity  of 
the  processed  image,  we  show  in  Figure  38  (a)  and  (b),  the  direct 
input  image  and  the  processed  output  image.  The  contrast  of  the  in¬ 
put  image  was  only  1%,  just  below  that  required  for  human  recogni¬ 
tion.  The  processed  image  has  a  much  improved  contrast  of  19*.  it 
is  clearly  visible  by  a  human  observer. 
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One  drawback  of  using  the  649F  plate  is  the  low  diffraction  ef-  ' 

ficiency  for  inputs  with  low  modulation.  When  the  intensity  of  the  I 

processed  image  is  too  low,  other  sources  of  noise  such  as  detector  • 

noise  and  ambient  noise  become  significant,  resulting  in  a  lower  out-  ; 

put  image  contrast.  A  common  technique  in  improving  the  efficiency 
of  recordings  made  with  photographic  materials  is  to  bleach  the  ; 

1 8  i 

emulsion  .  The  bleaching  operation  converts  the  developed  opaque  j 

silver  particles  into  transparent  silver  salts.  While  the  use  of  j 

bleaching  techniques  can  improve  the  diffraction  efficiency  of  the 
recording,  the  resulting  SNR  is  invariably  worse.  Thus,  bleaching 
would  not  be  desirable  in  our  present  application. 

The  use  of  a  material  such  as  dichromated  gelatin  can  signifi¬ 
cantly  increase  both  the  diffraction  efficiency  and  SNR  of  the  sys-  ; 

tern.  However,  dichromated  gelatin  is  far  too  insensitive,  particu¬ 
larly  to  red  light,  to  be  used.  Recently,  researchers  in  ERIM  de- 

19 

veloped  a  new  material,  the  silver  halide  gelatin  ,  which  exhibits 
a  diffraction  efficiency  and  noise  characteristic  almost  equal  to 
that  of  dichromated  gelatin.  The  sensitivity  of  the  material  how¬ 
ever,  is  comparable  to  that  of  a  photographic  plate  (optimum  exposure 
is  approximatel y  250  yJ/cm  ).  The  processing  procedure  for  the  ■ 

material  for  low  modulation  input  signals  has  not  yet  been  optimized. 

Nevertheless,  our  preliminary  experiments  indicate  that  a  substantial 

improvement  of  the  system  performance  can  be  achieved  by  substituting  j 

the  649F  plate  with  this  new  material.  In  Figure  37,  we  show  the 

experimental  results  obtained  with  the  silver  halide  gelatin  material  ; 

together  with  the  extrapolated  performance  curve.  We  see  that  the  ;■ 

system  performance  is  improved  and  the  processed  image  would  be  rec-  \ 

i 

ognizable  by  a  human  observer  if  the  input  has  an  image  contrast  of  j 

0.17%  or  higher.  A  positive  gain  in  image  contrast  can  be  obtained  j 

for  inputs  as  low  as  0.01%.  Moreover,  the  processed  image  obtained  j 

with  the  silver  halide  gelatin  material  is  significantly  brighter. 

For  example,  for  the  case  where  the  input  image  contrast  was  0.6%, 
the  intensity  of  the  processed  image  was  an  order  of  magnitude  higher 
than  that  obtained  with  the  649F  plates 
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5.  CONCLUSIONS  AND  RECOMMENDATIONS 

We  have  introduced  a  grating  interferometric  imaging  system  for 
imaging  through  scattering  media.  The  system  has  been  demonstrated 
to  be  capable  of  substantially  improving  the  image  contrast.  The 
system  is  shown  to  be  relatively  insensitive  to  normal  vibration  and 
air  disturbance,  making  it  feasible  for  field  applications. 

A  better  indication  of  the  system  capability  is  to  express  the 
performance  in  terms  of  visibility  range.  That  is,  the  maximum 
range  within  which  the  target  image  is  visible  through  a  scatter¬ 
ing  medium.  With  the  2%  minimum  contrast  criterion,  the  visibility 
range  is  related  to  the  visibility  attenuation  coefficient  ov  by 
the  relationship 

Rv  .  -  tn(o.O 2/.v) 

In  Figure  39,  we  plotted  the  visibility  range  curves  for  the  cases 
with  and  without  the  use  of  the  interferometric  system.  We  can  see 
that  the  visibility  range  is  increased  about  50%.  As  we  have  pointed 
out  before,  it  requires  more  than  a  2%  image  contast  to  recognize  a 
complex  target  (e.g.,  a  tank  within  a  cluster  of  rocks).  Let  us  as¬ 
sume  that  it  requires  a  5%  image  contrast  to  achieve  recognition  for 
a  particular  scene.  The  corresponding  visibility  ranges  are  also 
plotted  in  Figure  39.  We  find  that  with  the  5%  minimum  contrast 
criterion,  the  use  of  the  interferometric  system  would  improve  the 
visibility  range  by  80%.  That  is,  the  user  would  be  more  able  to 
detect  and  recognize  a  hostile  target  at  a  range  of  at  least  1.8  km 
compared  to  only  a  1  km  capability  otherwise.  In  general,  the  more 
complicated  the  target  scene,  the  larger  will  be  the  expected  im¬ 
provement  in  visibility  range  provided  by  grating  interferometer 
imaging  concept . 

The  use  of  an  imaging  grating  interferometer  for  imaging  through 
scattering  media  is  an  attractive  and  feasible  concept.  Supporting 
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FIGURE  39  VISIBILITY  RANGE  AS  A  FUNCTION  OF  VISIBILITY  ATTENUATION 
COEFFICIENT  OF  THE  SCATTERING  MEDIA 


imaging  experiments  performed  during  this  study  utilized  photographic 
materials  as  detectors.  For  near  real  time  operation,  a  sensitive 
real  time  recording  material  would  be  required.  One  possible  ap¬ 
proach  is  to  construct  an  image  converter-intensif ier  in  front  of  a 
real  time  light  modulator.  The  image  converter-intensif ier  would 
actively  increase  the  system  sensitivity  and  allow  the  use  of  IR  il¬ 
lumination.  An  alternate  approach  is  to  detect  the  stationary  in¬ 
tensity  pattern  produced  by  the  unscattered  light  with  a  CCD  photo¬ 
detector  array.  The  processing  could  be  performed  by  electronic 
(analog  or  digital)  filtering  and  subtraction  techniques.  Further 
work  is  needed  to  develop  the  concept  introduced  in  this  report  into 
a  practical,  field  usable  real  time  imaging  system  design. 
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